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ABSTRACT 

We present a study of the light curve of the transiting exoplanet CoRoT-2b, aimed at detecting the secondary eclipse and measuring its 
depth. The data were obtained with the CoRoT satellite during its first run of more than 140 days. After filtering the low frequencies 
with a pre-whitening technique, we detect a 0.0060+0.0020% secondary eclipse centered on the orbital phase 0.494±0.006. Assuming 
a black-body emission of the planet, we estimate a surface brightness temperature of Tp coRoT= 1910!^j;o K. We provide the planet's 
equilibrium temperature and re-distribution factors as a function of the unknown amount of reflected light. The upper limit for the 
geometric albedo is 0.12. The detected secondary is the shallowest ever found. 
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1. Introduction 

Out of the circa 350 exoplanets found to date, the 50 that tran- 
sit their host stars allow several follow-up studies, based either 
on the occultation of part of the stellar disk by the planet (tran- 
sits) or on the disappearance of the planet behind the star (sec- 
ondary eclipses). Observations of the secondary eclipses of sev- 
eral exoplanets have enabled measurements of their emitted flux 
at different wavel engths, pr oviding clues to t heir atmospheres 
jCharbonneau et al . 2005; Deming et alj|2005l) . Temperature in- 
versions at stratospheric levels have been found o n several occa- 
sions among the hottest known tra nsiting Jupiters dBurrows et alj 
I200 7'-, Knuts on et 311 12008', 2009), and several detections of the 
modulatio n of the planetary flux wi th orbital phase have been 
claimed bv lKnutson et al.1 (l2007l) and lSnellen et alj (l2009l) 

The transiting planet CoRoT-2b dAlonso et all 1200 8*) was 
disco vered during the fi rst ~150-d pointing of the CoRoT mis- 
sion dBaglin et al.L 12006 ). The spectral type of the star (GOV) 
and the period of the planet (1.7 d) su ggest that it be l ongs t o the 
pM class of exoplanets as defined by iFortnev et al.l (l2008l) . for 
which the strong incident flux is thought to cause stratospheric 
thermal inversions. 

One of the intriguing characteristics of the hot Jupiters is 
their variety in planetary radius. While for several planets the 
measured radius agrees with theoretical models of planetary for- 
mation and evolution, some of them appear larger than expected. 
A range of scenarios have been suggested to account for the 
phenomenon, but no clear picture has been delineated to date. 
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CoRoT-2b not only belongs to this class of "bloated" planets, 
but it is one with the most difficult radius as to explain, because 
of its high mass of S.lMjup. 

The host star of CoRoT-2b exhibits several signs of youth 
dBouchv et all l2008h . namely the detection of the Li I 6708 A 
absorption Une, an inversion of the cores of the Call H and 
K lines, and a fas ter-than-expected rotation period of 4.5 d 
(La nza et all l2009t) . While these "patterns" are traditionally at- 
tributed to stellar youth, it is also possible that the stellar evo- 
lution has been affected by tidal interactions with the close-in 
massive planet. Recentlv. Ijackson et all d2009h and lPontI d2009l) 
have suggested some observational evidence for strong tidal ef- 
fects on exoplanets' host stars. 

The activity of the star leaves footprints in the photometric 
light curve of CoRoT-2 that can be used to infer characteristics 
of the distribution and lifetimes of the spots on the stellar sur- 
face. By modeling the flux of the star in the parts of the light 
curve without transits, Lanza et al. I (l2009h observed a cyclic os- 
cillation of the total spotted area of the star with a period of 
28.9+4.3 days. In a different studv. lValio et all d2009 ) looked at 
the effects of the occultation of spots in the planet's path along 
the stellar surface. 

In this paper we describe detection of the secondary eclipse 
in the white light curve of the CoRoT public data. The tech- 
nique is simil ar to the one used to d etect the secondary eclipse 
of CoRoT-lb dAlonso et all l2009bh . A tentative detection of a 
5.5 x10'^ echpse, ba sed on the same data, has been claimed by 
Alon so et al.l d2009ab . Here we refine the analysis by carrying 
out a more careful filtering of the stellar variability, estimate the 
significance of the detection, and discuss the results. 
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2. Observations 

CoRoT-2b was observed during the first Long Run pointing 
of the satellite, which lasted 142-days. We used the data cor- 
rected to the N2 leve l (the processing steps are described in 
lAuvergne et alj l2009l) . which contains 369695 flux measure- 
ments, with a time sampling of 512-s from the first 5.2 days of 
data, after which it was changed to 32-s. The standard devia- 
tion of the normalized data after filtering the low frequencies as 
described below is 0.0014 (in units of normalized flux) for the 
32-s sampled data, while it is 0.00056 for the 512 s data. For 
comparison, the photon noise level in the 32-s data is of 0.001 1, 
revealing that the real data is only about 20% above the photon 
noise level, increasing to 90% above the photon noise in the 512- 
s part. An inspection of the 2MASS plates reveals that the star 
has a close companion 3.5 mags fainter in V, whose flux falls 
completely inside the aperture mask. We therefore subtracted a 
constant value of 5. 6% of the median le vel of the curve before 
normalizing it, as in lAlonso et al.l (l2008l) . 



3. Analysis of tlie CoRoT light curve 

The two most remarkable features of the light curve of CoRoT-2 
are the transits by its companion and the modulation caused by 
the stellar activity and rotation. Both features tend to hide any 
secondary eclipse in the light curve; therefore, in this section we 
explain our efforts to filter these two features, while avoiding to 
dilute or erase the signal of the secondary eclipse. 

To filter out the primary transit we could ignore the data 
points inside the transit. However, inclusion of data gaps at the 
orbital period of the transiting planet might result in spurious 
signals at this period and its harmonics, risking an effect on the 
signal of secondary eclipse. To reduce this risk, we subtracted 
the best fit solution of Alonso et al. (2008) from the curve, in- 
stead of cutting the parts of the light curves where the transits 
are detected. We plot the light curve with the tra nsits removed 
this way in Fig. [1] As noted bv lValio et al.l (l2009h . the residuals 
in the individual transits show the effect of occulted spots during 
the transit, so they might add some high frequencies in a very 
short phase of the planet's orbit. 

To filter out the stellar activity signal, we pr e-whitened the 
light curve using period04 (Lenz & Breger, l2005h . a well-known 
and tested technique for identifying periodic signals in variable 
stars. Basically, it consists of successively locating the highest 
amplitude in the Fourier power spectrum and fitting a combina- 
tion of sinusoids to the data that contain the frequency with the 
highest amplitude and all the previously identified frequencies. 
Before performing the period04 pre-whitening of the data, we 
averaged them into 30-min bins. We searched for frequencies 
between and 24 c/d (cycles per day). 

The amplitude spectra before and after filtering the low fre- 
quencies are plotted in Fig. |2] and the mean noise level in the 
0-24 c/d range in the filtered spectum is of 6.1x10"'', with a dis- 
persion of 4.2x10"*'. 

The search for the secondary eclip se was performed us - 
ing the same techniques as described in lAlonso et al.l (l2009bl) . 
Briefly, the search consisted of: 

- Selecting an orbital phase 0, where the secondary depth is 
going to be evaluated. 

- Performing low order polynomial fits (linear fits in this par- 
ticular case) to the parts of the light curve surrounding the 
phase (pi. The regions where the fits are evaluated are at a 
distance of slightly more than the transit duration, to avoid 
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Fig. 1. Light curve of CoRoT-2 where the modeled transits have 
been subtracted. For display purposes, the data have been binned 
in 30 min bins. 
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Fig. 2. The amplitude spectra of the CoRoT-2 light curve, be- 
fore (top) and after (bottom) pre-whitening of the 123 frequen- 
cies needed to model the low frequency content of the signal. To 
avoid including gaps in the data, t he transits \yere prev iously cor- 
rected by the best-fit model of Alonso et al.l ( l2009al) . The inset 
shows a zoom around the planet's orbital frequencies. The ver- 
tical solid and dashed green lines mark the planet's period and 
twice that value, respectively. 



including an hypothetical secondary centered on (/i,- in the 
fits. The total number of fits will thus be the same as the total 
number of observed secondary eclipses. 
- Fitting a trapezoid with the same duration and shape as ob- 
tained from the transit parameters, and centered on 0,-. The 
only free parameter is the trapezoid depth. 

The resulting plot of the significance and depth of a sec- 
ondary eclipse as a function of the orbital phase is given in Fig. [3] 
The highest peak in this diagram appears close to phase 0.5, cor- 
responding to a depth of about 6x10"^. 

To make sure that any of the corrections explained above 
did not add or remove some signal of the secondary eclipse. 
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Fig. 3. Depth of a trapezoid with the shape and duration of a 
transit, as a function of the planet's orbital phase. The maximum 
is centered on phase 0.5, coiTesponding to the secondary eclipse. 
In gray, the result in a light curve where a secondary eclipse of 
1x10"^ has been added to the raw light curve. 



we inserted a 0.01% secondary eclipse in phase 0.5 to the raw 
light curve, and performed the same analysis on this simulated 
data. The whole procedure, using the same filtering and repeat- 
ing the pre-whitening of the curve, was performed on the simu- 
lated data. The results are plotted in Fig.|3] along with the results 
of the real data. The highest peak in this case is also centered at 
0.5, and the fitted depth is of about 1.6x10"^, i.e., the sum of the 
inserted eclipse and the signal present in the data. We thus sug- 
gest that the whole preparation of the curve before performing 
the secondary search does not significantly alter the signal we 
want to measure. 

The phase-folded light curve during the phases of expected 
secondary eclipse is plotted in Fig. |4] together with the best fit- 
ted trapezoids. We explored the goodness-of-fit through func- 
tion around phase 0.5, in a grid of secondary eclipse centers 
and depths. The model used for the secondary was a trapezoid 
with duration and shape fixed to those of the transit. The ob- 
served secondary eclipse was binned into bins of 0.001 in phase 
(~2.5 min), and we estimated the error of each bin as the stan- 
dard deviation of the measurements inside the bin divided by the 
square root of the number of points inside the bin. We plot the 
resulting dependence and the formal sigma contours in Fig.|5] 
A secondary eclipse with a depth of 0.0060±0.0020% is detected 
in the CoRoT light curve at an orbital phase of 0.494+0.006. 

To estimate the significance of the secondary eclipse detec- 
tion, we performed the following test, which preserved the pres- 
ence of correlated noise in the data. From the data we subtracted 
the fitted secondary eclipse and shifted the residuals circularly, 
and then fitted a trapezoid with the phase and shape expected 
from the primary eclipse. Again, the only parameter was the 
secondary eclipse depth. All our 1000 trials resulted in depths 
less than 0.006%, the depth derived from the real data. The 
mean fitted depth was 0.0001+0.0008%, consistent with zero. 
Additionally, we checked the uncertainties in the depth measure- 
ment by two different methods. In the first one, we reinserted 
the fitted secondary eclipse after the residual curves of the pre- 
vious test were constructed and evaluated the fitted depths. We 
estimated the 1-cr uncertainty from the standard deviation of the 
depths that were fitted during this test. The result of this test gave 
a depth of 0.0059+0.0008%. As a second method, we fitted two 
Gaussian functions to (1) the distributions of points inside total 
secondary eclipse and (2) a subset of the part outside secondary 
eclipse, with the same number of points as (1) and with a ran- 
domly selected starting point. In the construction of the distribu- 
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Fig. 4. Phase-folded curve of CoRoT-2b during the phases of 
secondary eclipse. The data have been binned in 0.01 in phase 
(~25 min), and the 1-cr error bars in each bin calculated as the 
standard deviation of the data points inside the bin divided by the 
square root of the number of points inside the bin. The solid line 
shows our best-fit trapezoid, with a depth of 0.0060+0.0020%, 
and centered on phase 0.494+0.006. The dashed line shows a 
secondary eclipse centered on phase 0.5. 
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Fig. 5. The x'- value for different centers and depths of the sec- 
ondary eclipse, and the 1,2,3, and 4-cr confidence limits. 



tions, to minimize the eff'ect of an arbitrarily selected histogram 
bin size, we used random values of the sizes of the histogram 
bin, between 1x10"^ and 7x10"^ in units of normalized flux. 
We repeated this test using 500 diff'erent starting points and val- 
ues of the bin, and we estimated the 1-cr errors as above. This 
test resulted in a depth of 0.0066+0.0020%. We adopted a final 
value of the secondary eclipse depth of 0.006+0.002% from the 
results of the different tests. 

The secondary eclipse appears slightly off'set, at a 1-cr 
level, from the expected center for a circular orbit, which 
could be interpreted as a possible non-zero eccentricity of 
e cos <y=0.0 1+0.01. This eccentricity is within the 1-cr uncer- 
tainty of the orbital solutio n obtained from the radial velocity 
measurements (0.03+0.03 in lAIonso et"ani2008h . 
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4. Discussion 

In the optical, the observed flux of the planet can be expressed: 

— ^p,reflected ^p,i"eemitted "t" ^pjnternah (1) 

i.e., the sum of a component of reflected light, a thermal re- 
emission of the incident stellar flux, and a thermal emission not 
related to the incident flux from the star (e.g., internal heat or 
emission due to tidal forces). The reflected light can be esti- 

(R \^ 
where Ag is the geometric albedo. 

We can calculate the thermal component of the planet's emis- 
sion (the combin ation of f p.reemitted an d Fp.intemai) using the same 
procedure as in i Alonso et al.l (l2009bl) . For that purp ose, we used 
a value of Teff=5625+120 K for the star (Alonso eTall l2008l) 
and the model spectrum of a G8 V star from Pickles ( 19981) . We 
calibrated the model spectrum to obtain the same integrated flux 
as a Planck function with the Ten of the star. Taking the CoRoT 
spectral response function into account (" Auvergne et al.L l2009l) . 
the ratio of areas of the star and the planet, and the incident flux 
lost in the reflection (Ab/^* where Ab is the Bond albedo), we 
searched for the brightness temperature that matched the ob- 
served eclipse depth. We assumed a black-body emission for the 
planet and considered different values of the albedo. We assumed 
Ab - Ag, a reasonable assumption since more than 50% of the 
stellar flux is emitted at the CoRoT wavelengths. As an example, 
Ab ~ Ag within 20% for the four solar s ystem giant planets when 
integrating Ag between 0.4 and 0.9 ;um (lKarkoschkaill994 . The 
different possible solutions are plotted in Figure |6] The zero- 
albedo brightness temperature in the CoRoT bandpass calculated 
this way and including the uncertainties on the Tefr resulted in 
Tp,coRoT=1910;^^QQ K. If we further assume that internal - and 
that the planet is in thermal equilibriurrQ, this temperature favors 
high values of the re-distribution factor /=0. 60+0. 14 that are 
greater than the maximum expected for no re-distribution from 
the day to the night sides (/ - 0.5). This may be explained by (i) 
departures from a blackbody planetary emission, (ii) a non-zero 
albedo, or (iii) emission of significant internal energy from the 
planet. We defer the first possibility to future studies of radiative 
transfer in this planet. 

Assuming blackbody radiation and negligible Fp.intemab 
Fig. |6] shows that likely values of the albedo (corresponding 
to / between 0.25 and 0.5) ai-e Ag = 0.06 + 0.06, which con- 
firms the very low reflection of the planet obtained theoreti- 



2100 r 



cally(e.g.lMarlev et al.|[T999t ISudarskv et aT]|2000t ISeager et all 



I2OOOI: iHood et al]|2008l) . and observationally: upper limits for 
the albedos of exoplanets in similar conditions as CoRoT-2b 
have been reported by Charbonneau et al. (1999, Ag <0.3 for 
T Bootis b), Leight et al. ( 2003a,b, A^ <0.39 for t Bootis b, 
Ag <0.12 for HD 75289), bv lRoweet aLl(l2008l) on HD 209458b 
(a 3-cr upper limit of 0.17). Finally, CoRoT-lb was found to be 
such t hat Ag < 0.2 using CoRoT's red channe l of the light 
curve ([Snellen et all 12009) . and independently by lAIonso et al.l 
(I2009bl) ' using the white channel. 

To distinguish among the different components of the ob- 
served planetary flux, we would need chromatic information. 
The CoRoT data are delivered in three band-passes, but we could 
not reach the level of precision needed to detect the secondary in 
these data. We attribute this to the difficulty of achieving a good 
correction of the jitter in the satellite's pointing. 

As a concluding remark, the depth of the secondary eclipse 
of CoRoT-2b is roughly the same as a transit of a l/?e planet in 
front of a solar radius star. The significance of the detection of 
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Fig. 6. The equilibrium temperature of CoRoT-2b as a func- 
tion of the albedo matching the observed secondary eclipse 
depth. The zero-albedo equilibrium temperature is about 1 9 1 K. 
The dotted lines plot the re-distribution factors /. Realistic re- 
distribution factors require the contribution of some reffected 
light (non-zero albedo) or the presence of an additional source 
of energy in the planet (assumed zero in this plot). 



the secondary eclipse thus emphasizes the excellent capabilities 
of the CoRoT mission at detecting planets with radii of only a 

few Rg, . 
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